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Abstract 

We study multi-dimensional normal approximations on the Poisson space by means 
of Malliavin calculus, Stein's method and probabilistic interpolations. Our results yield 
new multi-dimensional central limit theorems for multiple integrals with respect to Pois- 
son measures - thus significantly extending previous works by Peccati, Sole, Taqqu and 
Utzet. Several explicit examples (including in particular vectors of linear and non-linear 
functionals of Ornstein-Uhlenbeck Levy processes) are discussed in detail. 
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1 Introduction 

Let {Z,Z,fj,) be a measure space such that Z is a Borel space and ^ is a o"- finite non- 
atomic Borel measure. We set = {B G Z : /i(-B) < oo}. In what follows, we write 
N = {N{B) : B G Z^} to indicate a compensated Poisson measure on (Z, Z) with control 
^. In other words, iV is a collection of random variables defined on some probability space 
(fl,J-,¥), indexed by the elements of Z^ and such that: (i) for every B,C £ Zf^ such that 
B D C = 0, the random variables N{B) and N(C) are independent; (ii) for every B G Z^, 

N{B) ' = ' N{B) - n{B), where N{B) is a Poisson random variable with paremeter n{B). A 
random measure verifying property (i) is customarily called "completely random" or, equiva- 
lently, "independently scattered" (see e.g. |24j). 

Now fix d > 2, let F = C L'^{a{N),F) be a vector of square-integrable 

functionals of N, and let X = {Xi, . . . be a centered Gaussian vector. The aim of this 
paper is to develop several techniques, allowing to assess quantities of the type 

dniF,X) = sup\E[g{F)]-E[g{X)]\, (1) 
gen 

where ^ is a suitable class of real-valued test functions on M'^. As discussed below, our princi- 
pal aim is the derivation of explicit upper bounds in multi-dimensional Central limit theorems 
(CLTs) involving vectors of general functionals of N. Our techniques rely on a powerful com- 
bination of Malliavin calculus (in a form close to Nualart and Vives |14|). Stein's method for 
multivariate normal approximations (see e.g. [H \T0\ [22] and the references therein), as well 
as some interpolation techniques reminiscent of Talagrand's "smart path method" (see [25], 
and also [3l[9]). As such, our findings can be seen as substantial extensions of the results and 
techniques developed e.g. in |51 [TUl [T7], where Stein's method for normal approximation is 
successfully combined with infinite-dimensional stochastic analytic procedures (in particular, 
with infinite-dimensional integration by parts formulae). 



The main findings of the present paper are the following: 



(I) We shall use both Stein's method and interpolation procedures in order to obtain explicit 
upper bounds for distances such as ([1]). Our bounds will involve Malliavin derivatives and 
infinite-dimensional Ornstein-Uhlenbeck operators. A careful use of interpolation techniques 
also allows to consider Gaussian vectors with a non-positive definite covariance matrix. As 
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seen below, our estimates are the exact Poisson counterpart of the bounds deduced in a Gaus- 
sian framework in Nourdin, Peccati and Reveihac |10j and Nourdin, Peccati and Reinert [U]. 

(II) The results at point (I) are applied in order to derive explicit sufficient conditions for 
multivariate CLTs involving vectors of multiple Wiener- Ito integrals with respect to N. These 
results extend to arbitrary orders of integration and arbitrary dimensions the CLTs deduced by 
Peccati and Taqqu [18j in the case of single and double Poisson integrals (note that the tech- 
niques developed in [18j are based on decoupling). Moreover, our findings partially generalize 
to a Poisson framework the main result by Peccati and Tudor |19) . where it is proved that, on 
a Gaussian Wiener chaos (and under adequate conditions), componentwise convergence to a 
Gaussian vector is always equivalent to joint convergence. (See also [10].) As demonstrated 
in Section 6, this property is particularly useful for applications. 

The rest of the paper is organized as follows. In Section 2 we discuss some preliminaries, 
including basic notions of stochastic analysis on the Poisson space and Stein's method for 
multi-dimensional normal approximations. In Section 3, we use Malliavin-Stein techniques to 
deduce explicit upper bounds for the Gaussian approximation of a vector of functionals of a 
Poisson measure. In Section 4, we use an interpolation method (close to the one developed 
in |9j) to deduce some variants of the inequalities of Section 3. Section 5 is devoted to 
CLTs for vectors of multiple Wiener-Ito integrals. Section 6 focuses on examples, involving 
in particular functionals of Ornstein-Uhlenbeck Levy processes. An Appendix (Section 7) 
provides the precise definitions and main properties of the Malliavin operators that are used 
throughout the paper. 

2 Preliminaries 
2.1 Poisson measures 

As in the previous section, (Z, Z, /x) is a Borel measure space, and iV is a Poisson measure on 
Z with control ^u. 

Remark 2.1 Due to the assumptions on the space {Z,Z,fi), we can always set (fi, J-", P) and 
N to be such that 



where 6z denotes the Dirac mass at z, and N is the compensated canonical mapping 



(see e.g. |20j for more details). For the rest of the paper, we assume that Q and A'^ have this 
form. Moreover, the cj-field is supposed to be the P-completion of the u-field generated by 
N. 

Throughout the paper, the symbol Lp'{^i) is shorthand for L^{Z, Z, fi). For n > 2, we write 
L?'{ijP') and L'^{ij,'^), respectively, to indicate the space of real- valued functions on which 
are square-integrable with respect to the product measure /i", and the subspace of L^(/i"') 




Lo{B)-fi{B), B£Z^, 00 gQ, 
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composed of symmetric functions. Also, we adopt the convention = Llip) = L'^{fJ,^) = 

and use the fohowing standard notation: for every n > 1 and every f,g £ L'^{fi^), 

f 1/2 
{f,9)L^ti") = f{zi,...,Zn)g{zi,...,Zn)^J.''{dzi,...,dZn), ||/||L2(^n) = {f,f)^2^^ny 

For every / G LP'{ijJ^), we denote by / the canonical symmetrization of /, that is, 

1 

/(Xi, . . . ,X„) = —^ l^f{Xa{l),- ■ • ,a;f7(n)) 

where a runs over the n! permutations of the set {!,..., n}. Note that, e.g. by Jensen's 
inequality, 

< II/IIl2(m") (2) 

For every / G n > 1, and every fixed z E Z, we write /(z, •) to indicate the function 

defined on Z'""""^ given by (zi, . . . , Zn~i) i— ?• f{z, zi, . . . , Zn~i)- Accordingly, /(z, •) stands for 
the symmetrization of the function /(z, •) (in (n — 1) variables). Note that, if n = 1, then 
/(z, •) = /(z) is a constant. 

Definition 2.2 For every deterministic function h G L'^{fi), we write Ii{h) = N{h) = 
h{z)N{dz) to indicate the Wiener-Ito integral of h with respect to N . For every n > 2 
and every f £ L'^{jj,"'), we denote by In{f) the multiple Wiener-Ito integral, of order n, 
of f with respect to N. We also set Inif) = In{f), for every f G L^(/i"), and Iq{C) = C for 
every constant C . 

The reader is referred e.g. to Privault [21j for a complete discussion of multiple Wiener-Ito 
integrals and their properties (including the forthcoming Proposition 12.31 and Proposition 12. 4p 
- see also [HIM] . 

Proposition 2.3 The following properties hold for every n,m > 1, every f £ L^(/i") and 
every g G L^fi"^): 

1. E[/„(/)] = 0, 

2. E[/„(/)/m(fii)] = n!(/,5-)i2(^n)l(„=^) (isometric property). 

The Hilbert space composed of the random variables with the form In{f), where n > 1 and 
/ G L1{ijP'), is called the nth Wiener chaos associated with the Poisson measure N. The 
following well-known chaotic representation property is essential in this paper. 

Proposition 2.4 (Chaotic decomposition) Every random variable F G L^(-F, P) = L^(P) 
admits a (unique) chaotic decomposition of the type 

oo 

F = E[F]+Y,In{fn) (3) 

n>l 

where the series converges in L^(P) and, for each n > 1, the kernel /„ is an element of Ll^fi^). 
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2.2 Malliavin operators 

For the rest of the paper, we shah use definitions and results related to Malliavin-type operators 
defined on the space of functionals of the Poisson measure N. Our formalism is analogous to 
the one introduced by Nualart and Vives [14]. In particular, we shall denote by D, 6, L and 
L~^, respectively, the Malliavin derivative, the divergence operator, the Ornstein-Uhlenbeck 
generator and its pseudo-inverse. The domains of D, 6 and L are written domD, domS and 
domL. The domain of is given by the subclass of L^(P) composed of centered random 
variables, denoted by Lg(P). 

Albeit these objects are fairly standard, for the convenience of the reader we have collected 
some crucial definitions and results in the Appendix (see Section [7|). Here, we just recall 
that, since the underlying probability space is assumed to be the collection of discrete 
measures described in Remark I2.H then one can meaningfully define the random variable 
U! 1— )• Fz{oj) = F{lo + 5z), to G fi, for every given random variable F and every z € Z, where 
5z is the Dirac mass at z. One can therefore prove that the following neat representation of 
D difference operator is in order. 

Lemma 2.5 For each F G domD, 

D^F = Fz- F, a.e.-fi{dz). 

A proof of Lemma [2.51 can be found e.g. in |14| [17] . Also, we will often need the forthcoming 
Lemma 12. 6| whose proof can be found in [T7] (it is a direct consequence of the definitions of 
the operators D, 5 and L). 

Lemma 2.6 One has that F G domL if and only if F & domD and DF G dom5, and in this 
case 

5DF = -LF. 

Remark 2.7 For every F G Lq(P), it holds that L~^F G domL, and consequently 

F = LL-^F = 5{-DL-^F) = -5{DL-^F). 

2.3 Products of stochastic integrals and star contractions 

In order to give a simple description of the multiplication formulae for multiple Poisson in- 
tegrals (see formula Q), we (formally) define a contraction kernel f -^l g on ^p+'J-''-' for 
functions / G L'^{fi^) and g G L'Kfj.'^), where p, g > 1, r = 1, . . . ,p A g and / = 1, . . . ,r, as 
follows: 

f*i aill,- ■ -^Ir-hh, ,tp_r,Sl, ,Sq-r) (4) 

fl'-{dzi,...,dzi)f{zi, , . . . ,Z/,7i,. . .,Jr-l,tl, , • • • ,tp-r) 

Xg{zi, , . . . ,Zi,Ji, . . . ,7^_/,Sl, , . . .,Sq-r). 

In other words, the star operator " -k^, " reduces the number of variables in the tensor product 
of / and g from p + q to p + q — r — I: this operation is realized by first identifying r variables 
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in / and (7, and then by integrating out / among them. To deal with the case / = for 
r = 0, . . . , p A q , we set 

f *r dill: • • • ) Irj tl: : ■ ■ ■ : tp—r: •Si, , • • • , Sq— r) 
= /(7I) • • • ,7r,il, , • • • ,tp-r)g{ll, • • • ,7r,Sl, , • ■ • ,Sq-r), 

and 

/*0 9{tl, ,...,tp,Sl,,...,Sq)= f 0g{tl, ,tp,Sl,,. . . ,Sq) = f{ti, , . . .,tp)g{si, , . . .,Sq). 

By using the Cauchy-Schwarz inequality, one sees immediately that f-k'^g is square-integrable 
for any choice of r = 0, . . . ,p A g , and every / S L^(/iP), g G 

As e.g. in [17| Theorem 4.2], we will sometimes need to work under some specific regularity 
assumptions for the kernels that are the object of our study. 

Definition 2.8 Let p > 2 and let f G L2(^p). 

1. The kernel f is said to satisfy Assumption A, if (/ ■kp~^ f) E L?'{ii^) for every r = 
l,...,p. Note that (fir'^f) e L'^ifiP) if and only if f e L^ijjP). 

2. The kernel f is said to satisfy Assumption B, if every contraction of the type 

{zi,...,Z2p-r-l) H> \ f\-kl \f\{zi,...,Z2p-r-l) 

is well-defined and finite for every r = l,...,p, every I = l,...,r and every {zi, Z2p-r-i) G 

^2p—r—l 

The following statement will be used in order to deduce the multivariate CLT stated in 
Theorem 15.71 The proof is left to the reader: it is a consequence of the Cauchy-Schwarz 
inequality and of the Fubini theorem (in particular, Assumption A is needed in order to 
implicitly apply a Fubini argument - see step (S4) in the proof of Theorem 4.2 in |17| for an 
analogous use of this assumption). 

Lemma 2.9 Fix integers p,q > 1, as well as kernels f G L1{^p) and g G Lg(/i'') satisfying 
Assumption A in Definition \2.8l Then, for any integers s,t satisfying l<s<t<pAg, one 
has that f -^l g ^ L?'{pP^'^~^~^), and moreover 

1. 

11/ 5'lli2(^p+9-t-s) = if *l-s f,g*'!i-sg)L^{fi^+'')^ 

( and, in particular, 

11/ /llL2(/i2p-s-i) = 11/ /||L2(^t+^) ); 

2. 

Il/*t 5111,2(^^+9-*-^) < ||/*p-* /||L2(^t+a) X fi(||^2(^t+«) 

= II/*? /llL2{/i2p-i,-t) X ll^*? 5'llL2{/i29-a-t). 
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Remark 2.10 1. Writing k = p + q — t — s, the requirement that 1 < s < t < pAq implies 
that 1^ — p| <k<p + q — 2. 

2. One should also note that, for every 1 < p < q and every r = 1, ...,p, 

[ {f4gfd,,P+^-^=f {f^i'-f){g*yg)di/, (5) 

for every / G L^^fi^) and every g S L^(/x''), not necessarily verifying Assumption A. 
Observe that the integral on the RHS of ([5]) is well-defined, since / •kp~^ / > and 
9 9>0. 

3. Fix p,q > 1, and assume again that / € L'^{fiP) and g £ L^(/i'') satisfy Assumption A in 
Definition 12.81 Then, a consequence of Lemma 12.91 is that, for every r = 0, ...,p Aq — 1 
and every / = 0, ...,r, the kernel f{z,-) g{z.,-) is an element of L^(/iP"'"''~*~^~^) for 
fj,{dz)-a\most every z G Z. 

To conclude the section, we present an important product formula for Poisson multiple 
integrals (see e.g. |6l |23] for a proof). 

Proposition 2.11 (Product formula) Let f G L^(/xP) and g € L^(/i''), p,q > 1, and sup- 
pose moreover that f -k^ g £ L^(/i^^''^''~') for every r = 1, . . . ,p A q and I = 1, . . . ,r such that 
I ^ r. Then, 

= E ( ^ ) ( ' ) E ( n lp^9) , (6) 

r=0 \ / V / ;=o ^ ^ 
with the tilde ~ indicating a symmetrization, that is, 

where a runs over all (p + q — r — l)\ permutations of the set {1, . . . ,p -\- q — r — 1} . 

2.4 Stein's method: measuring the distance between random vectors 

We write g G C'^(M'^) if the function : M*^ — )• M admits continuous partial derivatives up to 
the order k. 

Definition 2.12 1. The Hilbert-Schmidt inner product and the Hilbert - Schmidt 
norm on the class of d x d real matrices, denoted respectively by (•, ■)h.s. o,nd \\ ■ Wh.s.j 
are defined as follows: for every pair of matrices A and B, {A,B)h.s. '■= Tr{AB''") and 
\\-A\\h.s. = \/ (^1 ^)h.s.j where Tr{-) indicates the usual trace operator. 

2. The operator norm of a d x d real matrix A is given by ||^||op := sup|j^.|j^^^]^ ||^x||]Kd. 

3. For every function g : M.'^ ^ let 

,, ,, Igix) -g{y)\ 
\\g\\Lip ■= sup — II p — , 
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where \\ ■ ||]gd is the usual Euclidian norm on . If g ^ C"^(M'^), we also write 

\\Vg{x)-Vg{y)\y 



M2{g) := sup- 



F - y\ 



ifgeC' 



Msig) := sup 



\}lessg{x) — Kessg 



op 



F - y\ 



where Hess 5(2) stands for the Hessian matrix of g evaluated at a point z. 
4- For a positive integer k and a function g S C'^'(^'^) ; '^^ 

Qk 



b^'^'^lloo = max sup 
i<n<...<ifc<c(^g]]jd 



dxi-^ . . . dxi^ 



■9{x) 



In particular, by specializing this definition to g^"^^ = g" and g^^^ = g'" , we obtain 

52 



loo = max sup 

l<i\<i2<d 



\9 



b'"||oo = max sup 

l<h<i2<i:i<d -.fzwd 



dxi-^ dxi2 

Q3 



dxi-y^ ^3/^2 ^■^23 

Remark 2.13 1. The norm ||(7||Ljp is written Mi{g) in [1]. 

2. If 5 G Ci(R'^), then \\g\\Lip = sup \\Vg{x)\\js^d. If 5 G C2(M^), then 



-^2(5) = sup ||Hess5f(x) 



Definition 2.14 The distance d2 between the laws of two W^-valued random vectors X and 

Y such that E,\\X\\^d, E,\\Y\\j^d < 00, written d2{X,Y), is given by 

d2{X,Y) = sup\E[g{X)] -E[g{Y)]\, 
gen 

where Ti indicates the collection of all functions g S C^(M'^) such that \\g\\Lip ^ 1 cind M2{g) < 
1. 

Definition 2.15 The distance d^ between the laws of two W^-valued random vectors X and 

Y such that E||X||^j, E||y||^d < 00, written d3{X,Y), is given by 

d3{X,Y) = sup \E[g{X)] -E[g{Y)]\, 



where TL indicates the collection of all functions g G C'^(M'^) such that ||g"||oo < 1 o-nd 



a"'\\ < 

y Woo ^ 
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Remark 2.16 The distances d2 and ds are related, respectively, to the estimates of Section 
[3] and Section m Let j = 2,3. It is easily seen that, if dj{Fn,F) — )• 0, where Fn,F are 
random vectors in M'^, then necessarily Fn converges in distribution to F. It will also become 
clear later on that, in the definition of d2 and d^, the choice of the constant 1 as a bound 
for 1 1 5 1 1 Lip, M2{g), ||5'"||oo) ||s'"'||oo is arbitrary and immaterial for the derivation of our main 
results (indeed, we defined d2 and ^3 in order to obtain bounds as simple as possible). See the 
two tables in Section [4.21 for a list of available bounds involving more general test functions. 

The following result is a d-dimensional version of Stein's Lemma; analogous statements 
can be found in |51[Tni[22] - see also Barbour [Ij and Gotze [5j, in connection with the so-called 
"generator approach" to Stein's method. As anticipated, Stein's Lemma will be used to deduce 
an explicit bound on the distance d2 between the law of a vector of functionals of and the 
law of a Gaussian vector. To this end, we need the two estimates ([7]) (which is proved in |10]) 
and ([8]) (which is new). 

From now on, given a d x d nonnegative definite matrix C, we write ^^^(0, C) to indicate 
the law of a centered (i- dimensional Gaussian vector with covariance C. 

Lemma 2.17 (Stein's Lemma and estimates) Fix an integer d > 2 and let C = {C{i,j) : 
i,j = l,...,d}beadxd nonnegative definite symmetric real matrix. 

1. Let Y he a random variable with values in . Then Y ~ A/'(i(0, C) if and only if, 
for every twice differentiable function / : M'^ 1— t- M such that K\{C, Hess f{Y))H.s.\ + 
E|(y, V/(y))iR>d| < 00, it holds that 

E[(y, V/(y))K, - {C,Ressf{Y))H.s.] = 

2. Assume in addition that C is positive definite and consider a Gaussian random vec- 
tor X ~ A/'d(0,C). Let g : ^ ^ belong to C^(M'^) with first and second hounded 
derivatives. Then, the function UQ{g) defined by 

/•I 1 

Uog{x):= -E[g{Vtx + VT^tX)-g{X)]dt 

is a solution to the following partial differential equation (with unknown function f): 

g{x) - E[g{X)] = {x, Vf{x))^, - {C, Hess f{x))H.s. , x G M'^. 
Moreover, one has that 

sup \\RessUog{x)\\H.s. < \\C-^\\op\\C\\ll,^g\\up, (7) 

and 

MsiUog) < ^\\C-'\\%^ \\C\\opM2{g). (8) 

Proof. We shall only show relation ^ , as the proof of the remaining points in the statement 
can be found in [lOj. Since C is a positive definite matrix, there exists a non-singular symmetric 
matrix A such that = C, and A~^X ~ Md{0, Id)- Let Uog{x) = h{A~^x), where 

/■I 1 

h{x)= / -E[gA{Vix + VT^A-'X)-gA{A-^X)]dt 
Jo ^* 
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and gAix) = g{Ax). As A ~ Md{0,ld), the function h solves the Stein's equation 

{x, Vh{x))^, - Ah{x) = gA{x) - n9A{Y% 

where Y ~ A/'d(0, Id) and A is the Laplacian. On the one hand, as Hess gA{x) = A Hess g{Ax)A 
(recah that A is symmetric), we have 

M2igA) = sup ||Hess 5(A(a;) Hop = sup ||^Hess5(ylx)A||op 
= sup ||ylHess5(x)yl||op < \\A\\lpM2{g) 

= \\C\\apM2ig), 

where the inequahty above fohows from the well-known relation ||Ai?||op < ||^|| op ||i?||op. Now 
write hA-i{x) = h{A~^x): it is easily seen that 



HessC/o5(x) = Hess/iA-i(x) = A ^Ress h{A ^x)A ^ 



It follows that 



MsiUog) = MsihA-i) 

'|Hess /i^-i (x) — Hess /iA-i(yy Hop 



sup ■ 

xj^y \\x y\ 

I op 



sup ■ 



\A-^Ressh{A-^x)A-^ - A'^Ress h{A-^y)A-^^ 



xf^y \\x y 

^ ,, .-i,,2 \\Ress h{A~^ x) 

< L4 Mix sup -^^ ^ rf ^ 



, 2 ||Hess - Hess Hop || - 

"o „ X sup r. r. X — — . , . 

x^y W-y\\ \\A-^x-A-^y\ 



r, ||Hess/i(^ -"-x) — Hess ^yjwov _i 
< IIAMIlxsup — ^ —- X p Mlop 

xi^y 11^ X J\ y\ 



lop 



= \\C~^tl^^M^{K). 

Since Mj,{K) < ^^^M2{gA) (according to pi. Lemma 3]), relation ^ follows immediately. 



3 Upper bounds obtained by Malliavin- Stein methods 

We will now deduce one of the main findings of the present paper, namely Theorem 13.31 This 
result allows to estimate the distance between the law of a vector of Poisson functionals and 
the law of a Gaussian vector, by combining the multi-dimensional Stein's Lemma 12.171 with 
the algebra of the Malliavin operators. Note that, in this section, all Gaussian vectors are 
supposed to have a positive definite covariance matrix. 

We start by proving a technical lemma, which is a crucial element in most of our proofs. 

Lemma 3.1 Fix d> 1 and consider a vector of random variables F := {Fi, . . . , Fd) C -L^(P). 
Assume that, for all 1 < i < d, Fi ^ domD, and E[Fj] = 0. For all (p G C^(M^) with bounded 
derivatives, one has that 

d Q d 

D,^{Fu ...,Fd) = Y, Q^.<^iF){D,Fi) + R^AD^Fi, D,Fj), z E Z, 

i=l * *ii=l 



10 



where the mappings Rij satisfy 



\Rij[yi,y2)\ < sup 



t>{^)\ X \yiy2\ < :^\\(f>"\\oo\yiy2\- 

Proof. By the multivariate Taylor theorem and Lemma 12.51 

D,(P{F,,...,Fa) = ct){Fi,...,Fa){u: + 6,)-ct>{Fi,...,Fd){u;) 

= ^{Fi{oj + 6,),...,Fa{uj + 6,)) - <P{Fi{lo), . . .,Fa{oj)) 

= T5-0(i^iM, . . . ,Frf(a;))(Fi(a; + 5,) - F,{uj)) + R 
i=i 

d 



(9) 



i=l 



where the term R represents the residue: 



R = RiD.Fi, D,Fd) = R^ADzFi,D,Fj), 
and the mapping (2/1,^2) ^ Rij{yi,y2) verifies ©• ■ 

Remark 3.2 Lemma 13.11 is the Poisson counterpart of the multi-dimensional "chain rule" 
verified by the Malliavin derivative on a Gaussian space (see [8l[T0]). Notice that the term R 
does not appear in the Gaussian framework. 



The following result uses the two Lemmas 12 . 1 71 and 13 . H in order to compute explicit bounds 
on the distance between the laws of a vector of Poisson functionals and the law of a Gaussian 
vector. 

Theorem 3.3 (Malliavin-Stein inequalities on the Poisson space) Fix d > 2 and let 

C = {C{i,j) : i,j = 1, . . . ,d} be a d x d positive definite matrix. Suppose that X ~ Md{0, C) 
and that F = {Fi, . . . , Fd) is a -valued random vector such that E[Fj] = and Fi G douiD, 
i = 1, . . . ,d. Then, 

d2{F,X) < \\C~'\\op\\C\\ll,\ 



Yn{C{i,j) - {DF,,-DL-^F,)l2^^;, 



(10) 



^^\\C~'\\%'\\C\\op I Kdz)E 



\i=l ) \i=\ ) 



:iii 



Proof. If either one of the expectations in (jlOp and (|lip are infinite, there is nothing to prove: 
we shall therefore work under the assumption that both expressions (|10p -(|li p are finite. By 
the definition of the distance ^2, and by using an interpolation argument (identical to the 
one used at the beginning of the proof of Theorem 4 in |4j), we need only show the following 
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inequality: 



\E[giX)]-E[g{F)]\ 
<A\\C-'\UC\\y' 



(12) 



+l^B\\C-^%^C\\o, / Kdz)E 



^i=l / \i=l / 



for any g G C°°(M°') with first and second bounded derivatives, such that < A and 

M2{g) < B. To prove (|12p . we use Point (ii) in Lemma 12.171 to deduce that 



\E[g{X)]-E[g{F)]\ 

= \E[{C,BessUog{F))H.s. - {F,VUogiF))^,]\ 

i,j=l ■' k=l 

d 



C{i,j)j-^Uog{F) 

OXiOXj 



k=l 
d 



d 



6{DL-'Fk)—Uog{F) 



dxk 



k=l 



d 



D{—UogiF)y-DL-^F,^ 



d 

We write — — Uog{F) := (pkiFi, . . . ,Fj) = <f)}^{F). By using Lemma ISTTj we infer 



d 



DMFi,...,Fd) = Y,-Q^.^k{F){D,Fi) + R, 



1=1 



with Rk= Ri,j,k{DzFi,D:,Fj), and 



\Ri,j,k{yi^y2)\ < o sup 



02 



dxidxj 



■(pk{x] 



X \yiy2\- 
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It follows that 



\E[g{X)]-E[g{F)]\ 



dxidxj 



Uog{F) 



i,k=l 



92 



dxidxk 



iUo9{F)){DFi,-DL-^Fk)L2 



+ ^ E[{R,^,^k{DF,,DFJ),-DL-^Fk)L^^.)'_ 

i,j,k=l 



< JE[\\RessUog{F)\\ 



H.S.] X A E ^ [(^(^'•?') - {DF,,-DL~^F,)l2 



+ \R2\, 



where 



R2= Yl n{R^J,kiDF^,DF,),-DL-'Fk)L^^)]. 



i,j,k=l 



Note that ([7]) implies that ||Hess C/o5(-^)||//.5. < \\C ^\\op ||C'||op^||s'||Lip- By using 1^ and the 
fact 1 1 g'" 1 1 00 ^ M^^g), we have 



\Ri,j,k{yi^y2)\ < sup 



Q3 



dxidxjdxk 



Uoigiy)) 



X |yiy2| 



27r 



< _-M2(5)||C''|i;!/'||C|Up X \y,y2\ < -^B\\C-^%^C\\op x ly.y^l 

from which we deduce the desired conclusion. ■ 

Now recall that, for a random variable F = N{h) = Ii{h) in the first Wiener chaos of N, 
one has that DF = h and L~^F = —F. By virtue of Remark I2.16| we immediately deduce 
the following consequence of Theorem 13.31 



Corollary 3.4 For a fixed d >2, let X ^ A/rf(0, C), with C positive definite, and let 

F„ = {Fn,l, ...,Fn,d) = {N{hn,l), N {hn,d)) , U > 1, 

be a collection of d- dimensional random vectors living in the first Wiener chaos of N. Call Kn 
the covariance matrix of Fn, that is: Kn{i,j) = E[iV(/i„.j)A^(/i„.j)] = {hn,i,hnj) i2f^^y Then, 



d2{Fn,X) < ||C-i||,p||C7||i/2||C-ir„||^,.5.+ 
In particular, if 



dW2TT 



d „ 

\C'^\\%^C\\opY, / \hnAz)fKd^ 
i=l -^^ 



Kn{i,j) ^ C{i,j) and / |/i„,i(z)p/u((iz) ^ 



(13) 

(as n — 7- oo and for every i,j = l,...,d), then d2{Fn,X) — )• and Fn converges in distribution 
to X. 
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Remark 3.5 1. The conclusion of Corollary I3.4l is by no means trivial. Indeed, apart from 
the requirement on the asymptotic behavior of covariances, the statement of Corollary 
13.41 does not contain any assumption on the joint distribution of the components of the 
random vectors Fn- We will see in Section 5 that analogous results can be deduced for 
vectors of multiple integrals of arbitrary orders. We will also see in Corollary 14.31 that 
one can relax the assumption that C is positive definite. 

2. The inequality appearing in the statement of Corollary 13.41 should also be compared 
with the following result, proved in |10) . yielding a bound on the Wasserstein distance 
between the laws of two Gaussian vectors of dimension d > 2. Let Y ~ Md{0, K) and 
X ~ A/'(i(0,C), where K and C are two positive definite covariance matrices. Then, 
dw{Y,X) < Q{C,K) X \\C-K\\h.s., where 

QiCK) ■.= mm{\\C-'\\op \\C\\l/,M\K-'\\op \\K\\l/,^}, 

and dw denotes the Wasserstein distance between the laws of random variables with 
values in W^. 



4 Upper bounds obtained by interpolation methods 
4.1 Main estimates 

In this section, we deduce an alternate upper bound (similar to the ones proved in the previous 
section) by adopting an approach based on interpolations. We first prove a result involving 
Malliavin operators. 

Lemma 4.1 Fix d > 1. Consider d + 1 random variables Fi S L^(P), < i < d, such that 
Fi S domZ) and K[Fi\ = 0. For all g € C^(M'^) with bounded derivatives, 



ngiFi, 

where 



Fd)Fo]=E 



d 



Q^giFi, • • • , Fd){DFi, -DL-'Fo)l2 



.1=1 



+E [{R,-DL-^Fo)l2 



(m)J 



\E[{R,-DL-'Fo)L2i^ 



(14) 



1 

< -max sup 

2 id .rPlRd 



dxidxj 

Proof. By applying Lemma |3.1 



■gix) 



X / fM{dz)E 



Y,\DzFk\] \D,L-'Fo 



\k=l 



E[5(Fi,...,Frf)Fo] 
E[{LL-^Fo)g{Fi,...,Fd)] 
-E[S{DL-^Fo)g{Fi,...,Fd)] 
E[{Dg{Fi,...,Fd),-DL~^Fo)L2f^^,- 



E 



d 



^ Q^giFi, Fd){DF,, -DL-'Fo)l2 



1=1 



+ E[{R,-DL-'Fo)l2^^)] 



14 



and E[(i?, —DL ^i^o)L2(^)] verifies the inequality ([T] 

As anticipated, we will now use an interpolation technique inspired by the so-called "smart 
path method", which is sometimes used in the framework of approximation results for spin 
glasses (see [25]). Note that the computations developed below are very close to the ones used 
in the proof of Theorem 7.2 in [9]. 



Theorem 4.2 Fix d > 1 and let C = {C{i,j) : i,j = l,...,d} be a d x d covariance 
matrix (not necessarily positive definite). Suppose that X = [Xi, ...,X(i) ~ A/'rf(0, C) and 
that F = {Fi, . . . , Fd) is a W^-valued random vector such that E[Fj] = and Fi G dovnD, 
i = 1, . . . , d. Then, 



dz{F,X) < 



mCiiJ) - {DF,,-DL-^F,)l2^^)Y] 



(15) 
(16) 



Proof. We will work under the assumption that both expectations in (|15p and ()16p are finite. 
By the definition of distance d^, we need only to show the following inequality: 



1 



|E[<^(X)]-E[</)(F)]| < -||0"|U5^E[|C(i,j)-(^i^.,-^i"'^.-)L2 



/i(dz)E 



\i=l J \i=l J 



for any (j) G C (M ) with second and third bounded derivatives. Without loss of generality, we 
may assume that F and X are independent. For t E [0, 1], we set 



We have immediately 



^{t) = E[0(^/r^(Fi, . ..,Fd) + ViX)] 



1^(1) -^(0)1 < sup |^'(t)|. 

te{o,i) 



Indeed, due to the assumptions on 0, the function 1 1— )• ^{t) is differentiable on (0, 1), and one 
has also 



V(t) = ^E 



1 



d 

dxi 



Vl^t{Fi,...,Fd) + VtX 



:2l 



2Vt 2^1^ 



55. 
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On the one hand, we have 



2t 



d 



i=l 



E 



-^q^iVT^a + VtX)Xi 
dxi 



\a={Fi,...,Fa)^ 



E 



dxidxj 



- ta + VtX) 



dxidxj 



|a={Fi,...,Fd). 

■cl){VT^t{Fi,...,Fd) + ViX) 



On the other hand, 



« = ^E 

i=l 

= E^^ 



j=i 



E 



d 

dx. 



■(PiVT^t{F,,...,Fd) + Vtb)Fi 



\b=X 



We now write 4>^ (•) to indicate the function on defined by 

<Af (Fi, . . . , Frf) = —^{VT^t{Fi, ...,Fa) + Vtb) 

By using Lemma l4.ll we deduce that 
E[4\Fi,...,Fa)Fi] 

d 



E 



+ e[{ri-dl-'f,)l2 



where Rl is a residue verifying 

|E[(i?^,-DL-iF,)i2(^)]| 



< - max sup 



dxkdxi 



^\D,Fj\ I \D,L-'Fi\ 
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Thus, 

d 



E 



52 



dxidxj 



■0(^/r^(Fi, ...,Fd) + Vtb){DFi, -DL-^Fj 



\h=x 



+ 



52 



2=1 



dxidxj 



cPiVT^t{Fi, ...,Fd) + VtX){DFi, -DL-'F, 



.^E[E[(i?^,-Z)L-iF,)i2(^)] 



2=1 



|f)=X 



Putting the estimates on 21 and !B together, we infer 



1 



92 



dxidxj 



■cPiVr~t{Fu ...,Fd) + VtX){C{i,j) - {DFi, -DL-'Fj)l2^^)) 



We notice that 



1 

=J]e[e [{RI-DL-'F,)l2^^)] 



\b=X 



dxidxj 



■<l){Vl^t{Fi,...,Fd) + Vth) 



< 



and also 



52 



X 



53 



dxidxkdxi 



4{VT^{Fi,...,Fd) + Vtb) 



< (i-t)iinioo. 

To conclude, we can apply inequality ()17p and deduce the estimates 

|E[</.(X)]-E[,/.(F)]| 

< sup 

*e(o,i) 

1 

< -||(/>"|U ]E[|C(z,j) - {DF„-DL-'F,)l2^^)\] 



+ 



1-t 



fi{dz)E 



^\D,F,\] (Y,\D.L-'F,\ 



\i=l 



< - 



^ E[(C(i,i) - {DF„ -DL-^F^)l2(^)Y] 



fi{dz)E 



^i=l / \i=l / 



17 



thus concluding the proof. ■ 

The following statement is a direct consequence of Theorem I4.2[ as well as a natural 
generalization of Corollary 13.41 

Corollary 4.3 For a fixed d> 2, let X A/d(0, C), with C a generic covariance matrix. Let 
Fn = {Fn,i,...,Fn,d) = {N{K,i),...,N{K,d)), n>l, 

be a collection of d- dimensional random vectors in the first Wiener chaos of N , and denote 
by Kn the covariance matrix of Fn- Then, 

i=l •' ^ 

In particular, if relation M^) is verified for every i,j = 1, ...,d (as n ^ oo), then d^{Fn,X) — )• 
and Fn converges in distribution to X. 

4.2 Stein's method versus smart paths: two tables 

In the two tables below, we compare the estimations obtained by the Malliavin-Stein method 
with those deduced by interpolation techniques, both in a Gaussian and Poisson setting. Note 
that the test functions considered below have (partial) derivatives that are not necessarily 
bounded by 1 (as it is indeed the case in the definition of the distances ^2 and d^) so that the 
L°° norms of various derivatives appear in the estimates. In both tables, d > 2 is a given posi- 
tive integer. We write (G, Gi , . . . , Gd) to indicate a vector of centered Malliavin differentiable 
functionals of an isonormal Gaussian process over some separable real Hilbert space (see 
for definitions). We write {F,Fi, ...,Fd) to indicate a vector of centered functionals of N, 
each belonging to domZ). The symbols D and stand for the Malliavin derivative and the 
inverse of the Ornstein-Uhlenbeck generator: plainly, both are to be regarded as defined either 
on a Gaussian space or on a Poisson space, according to the framework. We also consider 
the following Gaussian random elements: X ~ A/'(0, 1), Xq ~ J^di^j C) and Xm ~ ^/^(O, M), 
where C is a, dx d positive definite covariance matrix and M is a d x d covariance matrix (not 
necessarily positive definite). 

In Table 1, we present all estimates on distances involving Malliavin differentiable ran- 
dom variables (in both cases of an underlying Gaussian and Poisson space), that have been 
obtained by means of Malliavin-Stein techniques. These results are taken from: |8j (Line 1), 
|10j (Line 2), |17j (Line 3) and Theorem 13.31 and its proof (Line 4). 

In Table 2, we list the parallel results obtained by interpolation methods. The bounds in- 
volving functionals of a Gaussian process come from |^ , whereas those for Poisson functionals 
are taken from Theorem 14.21 and its proof. 

Observe that: 

• in contrast to the Malliavin-Stein method, the covariance matrix M is not required to 
be positive definite when using the interpolation technique, 

• in general, the interpolation technique requires more regularity on test functions than 
the Malliavin-Stein method. 
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Table 1: Estimates proved by means of Malliavin-Stein techniques 



Regularity of 
the test function h 


Upper bound 


\\h\\up is finite 


|E[/?,(G)] - ^[h{X)]\ < 
\\hh,, ^E[{1 - {DG.-DL-HJ),, f] 


|/i is finite 


\E[hiGi,...,Gd)]-E[hiXc)]\ < 
||/i||L^p||C-i|U||C||i/yEt=iE[(C(i, j) - {DGi,-DL-^Gj)s^y] 


||/i||Lip is finite 


\E[h{F)]-E[h{X)]\ < 
\\h\\Lrpi^E[{l - {DF,-DL-^F)L2^^)y] 
+ J^I,{dz)E[{\D,F\r\D,L-'F\]) 


h G C2(M'^) 
|/i liip is finite 

M2{h) is finite 


\E[h{Fi,...,Fd)]-E[h{Xc)]\ < 
+M2{h)^\\C-m('\\G\\opJzKdz)E (^El^.i^il) (j:\D,L-^F,\^ 



"JVibk" 2: Estiniatcs i")rcn'cxi b\' means of inlerpolatioiis 



Regularity of 
the test function (f) 


Upper bound 


(j) G C2(R) 
10" loo is finite 


|E[0(G)]-E[0(X)]|< 
\\\ct>"U^E[{l - {DG-DL-^G)^Y\ 


10" loo is finite 


|E[0(Gi,...,Gd)]-E[0(XM)]| < 
W\U\jT.ij=inmhj) - {DGi,-DL-^Gj)s,n 


(j) G C3(R) 
10" loo is finite 
10"' loo is finite 


|E[0(F)]-E[0(X)]|< 
i||0"||ooy/E[(l - {DF-DL-^F)^.^^)^] 
Woo Jzf^idz)E[i\D,F\n\D,L-^F\)] 


G C^{W^) 
10" loo is finite 

10"' loo is finite 


|E[0(Fi,...,F,)]-E[0(Xm)]| < 


W\\oo^Ei,j=inmhj) - {DF„-DL-^Fj)L2^^)f 
+W'\\ooIzKdz)E (eJD^F.I^ (j:^\D,L-^F,\^ 
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5 CLTs for Poisson multiple integrals 



In this section, we study the Gaussian approximation of vectors of Poisson multiple stochastic 
integrals by an application of Theorem 13.31 and Theorem 14.21 To this end, we shall explicitly 
evaluate the quantities appearing in formulae pO|) - (jlip and (I15p -(|16 p . 

Remark 5.1 (Regularity conventions) From now on, every kernel / S L'^{fi^) is supposed 
to verify both Assumptions A and B of Definition 12.81 As before, given / G and for 

a fixed z ^ Z, we write f{z, •) to indicate the function defined on as (zi, . . . , Zp^i) i— )■ 

f{z,zi,... ,Zp-i). The following convention will be also in order: given a vector of kernels 
(/i) fd) such that fi S Ll{fi'P'), i = 1, d, we will implicitly set 

fi{z,-) = 0, z = 

for every z & Z belonging to the exceptional set (of ^ measure 0) such that 

U{z,-)4f,{z,.)iL\i,^^^P^~^''-^) 
for at least one pair (i, j) and some r = 0, ...,pi Apj — 1 and / = 0, r. See Point 3 of Remark 

Eini 

5.1 The operators Gf'^ and cf^ 

Fix integers p,q > and !(? — p| < k < p + q, consider two kernels / G Ll^fiP) and g G 
L1{^'^), and recall the multiplication formula ([6]). We will now introduce an operator G^''^, 
transforming the function /, of p variables, and the function g, of q variables, into a "hybrid" 
function G^''^{f,g), of k variables. More precisely, for p,q,k as above, we define the function 
{zi, . . . , Zk) I—)- G^''^{f,g){zi, . . . , Zk), from Z^ into M, as follows: 

Gl'\f,g){z,,...,z,) = X:E1(p+,-.-/=;.)H ( f ) ( r ) ( ^^^^ 

where the tilde ~ means symmetrization, and the star contractions are defined in formula 
and the subsequent discussion. Observe the following three special cases: (i) when p = q = 
k = 0, then / and g are both real constants, and Gq''^(/, gf) = f x g, (ii) when p = q > 1 and 
k = 0, then GQ'^{f,g) = p^f, g)L'^[^), (hi) when p = k = and q > (then, / is a constant), 
G^'^if, g){zi, Zq) = / X g{zi, Zg). By using this notation, ^ becomes 

p+q 

W)h{9)= E h{Gl\f,g)). (19) 

k=\q-p\ 

The advantage of representation (|19p (as opposed to ^) is that the RHS of (jl9p is an orthog- 
onal sum, a feature that will greatly simplify our forthcoming computations. 

For two functions / G L'^{fiP) and g € L^{fj,'^), we define the function {zi,...,Zk) ^ 
G^''^{f, g){zi, . . . ,Zk), from Z^ into M, as follows: 

Gf^(/,5)(-)= / Kdz)Gl-''''-\f{z,.),g{z,.)), 
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or, more precisely, 

G^'^if,g)izu...,zk) 



r=0 1=0 



E E ^)(* - ( * - 1 ) ( ' - 1 ) ( ' - 1 ) ■ ■ ■ ' ^'"^ 

Note that the impHcit use of a Fubini theorem in the equahty ()20p is justified by Assumption 
B - see again Point 3 of Remark 12.101 

The following technical lemma will be applied in the next subsection. 

Lemma 5.2 Consider three positive integers p, q, k such that p,q > 1 and |(? — p| V 1 < A; < 
p + q — 2. For any two kernels f € L^{fi^) and g G Ll{ijfl), both verifying Assumptions A and 
B, we have 



"'^'^ t=i 
where s{t, k) = p + q — k — t for t = 1, . . . ,p A q. Also, C is the constant given by 



(21) 



pAq 



(t-i) 



t=i 

Proof. We rewrite the sum in ()20p as 

pAq 



p-l \ ( q-l \ ( t-l 
t-1 J \t-l J \ s{t,k)-l 



-\ 2 



Cfc ^(/) 9)(^i, • • • , ^fc) = E ^t^i<s{t,k)<tf 9{zi, Zk), (22) 

t=i 



witha,=(*-i).( ^:;)(?:;)(,„^-/_i i, i<<<,-A,^Th„B. 



Zk 



df,\Gl'''{f,g)izu...,Zk)f 



< 



/pAq ^ \ ^ 

Jzk '^^^ \ Y^^^t'^'^<s{t,k)<tf -4^^'^'' 9{Z1, • • ■,Zk)\ 
/pAq \ /pAq \ 

(E«?) (E(ii<K*,^)<*/^?^''^5(-i'---'-'^))' 

pAq 

•"V / dp^\^<,i^^k)<t{f 4^'^^ 9{.z\,---,ZkS? 

pAq ■ 
C'Y'i-l<s{t,k)<t\\f *t ' 5lli2(^fc), 



t=l 
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with 

pAq pAq 

t=i t=i 

Note that the Cauchy-Schwarz inequahty 




has been used in the above deduction. ■ 
5.2 Some technical estimates 

As anticipated, in order to prove the multivariate CLTs of the forthcoming Section I5.3| we 
need to estabhsh exphcit bounds on the quantities appearing in (|1U|) - (|11|) and (|15p -(|16 p . in 
the special case of chaotic random variables. 



Definition 5.3 Let the integers p,q > 1 be such that max(p, q) > 1. The kernels f € L^{fjP), 
g G L^(/Li'^) are said to satisfy Assumption C if, for every k = \q — p\ V 1, . . . ,p + g — 2, 



{Gl~''''~\fiz,.),giz,.Wdf^'' 



fi{dz) < oo. 



(23) 



Remarlc 5.4 By using p8|) . one sees that (|23p is implied by the following stronger condition: 
for every k = \q — p\\/ 1, . . . ,p + q — 2, and every (r, /) satisfying p + q — 2 — r — I = k, one has 



Zk 



{f{z,-)4g{z,-Wdf,'^ 



n{dz) < oo. 



(24) 



One can easily write down sufficient conditions, on / and g, ensuring that ()24p is satisfied. For 
instance, in the examples of Section 6, we will use repeatedly the following fact: if both / and 
g verify Assumption A, and if their supports are contained in the finite union of rectangles of 
the type B x . . . x B, with fJ,{B) < oo, then is automatically satisfied. 

Proposition 5.5 Denote by the pseudo-inverse of the Ornstein- Uhlenbeck generator (see 
the Appendix in Section]^, and let F = Ip{f) and G = Iq{g) be such that the kernels f G 
L'^{ij,p) and g £ L'^^fj,'^) verify Assumptions A, B and C. If p ^ q, then 

n{a-{DF,-DL-'G)L2(^)f] 

P+<?-2 

< a'+p' d^^HG''^^f,g)f 

, I I Jzk 

k=\q-p\ 

p+q— 2 pAq 



< a^ + Cp'^ Y '^•Z]^ 



l<s(t,k)<t\ 



s{t,k) 



k=\q-p\ t=l 
^ P+q-2 pAq 

k=\q-p\ t=\ 
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If p = q, then 



E[ia-{DF,-DL-'G)L^^)r] 



k=i 

2p-2 pAq 



k=l t=l 



< 



a 



2p-2 pAq 



fc=l t=l 



where s{t,k) = p + q — k — t for t = 1,... ,p Aq. Finally, the constant C is given by 



pAq 



t=l 



(t-iy. 



p — 1 



t-l J \ t-l J \ s{t,k) - 1 



t-1 



Proof. We select two versions of the derivatives DzF = plp-i{f{z, •)) and DzG = qlq-i{g{z, •)), 
in such a way that the conventions pointed out in Remark 15. II are satisfied. By using the def- 
inition of and (1191). we have 



{DF,-DL-'G)l2(^^^ = {Dl,{f),q-'DI,{g))L2(^^) 



p p{dz)Ip^i{f{z,-))Iq-l{g{z,-)) 

Jz 

k=\q-p\ 



Notice that for i ^ j, the two random variables 



p{dz)h{Gl^^''-\f{zr),g{zr)) and / ^,{dz)I,{G^-^^'^-\f{z,•),g{z,■))) 



are orthogonal in L^(P). It follows that 

E[{a-{DF,-DL-^G)L2(^)) 

p+q~2 

= a^+p^ ^ E 

k=\q~p\ 

for p ^ q, and, for p = q, 

E[ia-{DF,-DL-'G)L2i^)f] 

2p-2 

= {pl{f,g)m^P)-af+p^Y.^ 

k=l 



(25) 



f^{dz)h{Gl-''''-\fiz,-),g{z,-))) 



(26) 



p{dz)h{Gt'''^\f{z,-),g{z,-))) 
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We shall now assess the expectations appearing on the RHS of (|25p and (|26p . To do this, fix 
an integer k and use the Cauchy-Schwartz inequality together with (|23p to deduce that 



p— i,<j— 1 / f( I 



< / ^,{dz) / /.(dz')VK[/i(Gr''^~'(/(^,-),5(^,-)))]VE[/i(G'r^'''~^(/(z',.),5(^',-)))] 



„{dz)\l I^J^HGl-'''-\fiz,-),g{z,W 



f,idz')J j^^ df,HGl-'''^\f{z', ■),giz', 



kl 



f^{dz)J I dt,'^{Gl-'''''\f{z,-),g{z,-W 



< oo. 



(27) 



Relation (|27p justifies the use of a Fubini theorem, and we can consequently infer that 



E 



f,{dz)h{Gl-'^''-\f{z,.),g(z,-))) 



f,{dz) / M(dz')E[4(Gr'''~'(/(^r),5(^,-)))4(Gr^'^-^(/(z',-),5(^',-)))] 
Jz 



k\ / ^i{dz) / ^i{dz') 



k\ / V 



zk 

~(P-l,g-l 



= k\ df,''iGl'''if,g)r. 

The remaining estimates in the statement follow (in order) from Lemma 15.21 and Lemma 12.9 
as well as from the fact that ||/||L2(^n) < ||/||L2(^n), for all n > 2. ■ 



The next statement will be used in the subsequent section. 



Proposition 5.6 Let F = (Fi, . . . , F^) := (/^^(/i), . . . , /^^(/d)) be a vector of Poisson Junc- 
tionals, such that the kernels fj verify Assumptions A and B. Then, by noting := minjgi, q^}, 



n{dz)E 



Ki=l J \i=l / 



Qi 6-1 



^ ;r E i^V^'?^ - ^ E E li<a+f'<2,.-i(a + h - iy}l\q, - a - 1)! 



; - 1 - a 



lL2(At«) 

Qi - 1- a 
Qi - b 



b=l a=0 

f '^h f\\L^^{^?1i-'^-*'A■ 
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Proof. One has that 



n{dz)E 



q* Jz 



< — [ n{dz)E 



To conclude, use the inequahty 

/ fi{dz)\E[DJ,{f)]f 
Jz 



q 6-1 



6=1 a=0 



- 1 

I- a 



q — 1 — a 



which is proved in [17, Theorem 4.2] (see in particular formulae (4.13) and (4.18) therein). 



5.3 Central limit theorems with contraction conditions 

We will now deduce the announced CLTs for sequences of vectors of the type 



i^^") = {Ft\. . . ) := (/,,(/r), . . .JqAfD), n > 1 



in). 



(28) 



As already discussed, our results should be compared with other central limit results for 
multiple stochastic integrals in a Gaussian or Poisson setting - see e.g. [8l [TOl IT2| 1131 [T8l 119). 
The following statement, which is a genuine multi-dimensional generalization of Theorem 5.1 
in |17| . is indeed one of the main achievements of the present article. 

Theorem 5.7 (CLT for chaotic vectors) Fix d>2, let X J\f{0, C), with 

C = {Cii,j):iJ = l,...,d} 

a d X d nonnegative definite matrix, and fix integers qi,...,qd > 1. For any n > 1 and 
i = 1, . . . , d, let f^"'^ belong to L^(^'^'). Define the sequence {F*^"^ : n > 1}, according to 
and suppose that 



hm E[i^(")Ff )] = X hm (/f^ /r)^^{/'«) = ^(^'•^■)' ^<hj<d. (29) 

Assume moreover that the following Conditions 1-4 hold for every k = l,...,d: 



An) AnX 
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1. For every n, the kernel fj^^ satisfies Assumptions A and B. 

2. For every I = l,...,d and every n, the kernels fj,"'^ and fj^'^^ satisfy Assumption C. 

3. For every r = 1, . . . , g^, and every / = 1, . . . , r A {q^ — 1), one has that 

II f (") ,/ f (")|| > n 

as n ^ OO. 



I Asn^ oo, f^,, dfii^ i^fl!'') ^ 0. 

Then, F^") converges to X in distribution as n ^ oo. The speed of convergence can be assessed 
by combining the estimates of Proposition 15.51 and Proposition 15.61 either with Theorem 
(when C is positive definite) or with Theorem^J^ (when C is merely nonnegative definite). 



Proof. By Theorem 14.2 



^E[iCiz,j) - {DFi''\-DL-^F^^^)L.^^))^] (30) 



d3iF(^\X)<^ 

+i/^,(^.)E {E\D.Ft'\)\E 



d.l-^f!-''^] 



(31) 



so that we need only show that, under the assumptions in the statement, both ()30p and (|3ip 
tend to as n — )• (X). That (|30p tends to is a direct consequence of the estimates in Propo- 
sition whereas Proposition 15.61 shows that (|3ip converges to 0. This concludes the proof. 

If the matrix C is positive definite, then one could alternatively use Theorem 13.31 instead 
of Theorem [ 



Remark 5.8 Apart from the asymptotic behavior of the covariances (j29p and the presence 
of Assumption C, the statement of Theorem 15.71 does not contain any requirements on the 
joint distribution of the components of F^^\ Besides the technical requirements in Condition 
1 and Condition 2, the joint convergence of the random vectors F^"^ only relies on the 'one- 
dimensional' Conditions 3 and 4, which are the same as condition (II) and (III) in the statement 
of Theorem 5.1 in [17j. See also Remark 13.51 

6 Examples 

In what follows, we provide several explicit applications of the main estimates proved in the 
paper. In particular: 

• Section 16.11 focuses on vectors of single and double integrals. 

• Section [6^ deals with three examples of continuous-time functionals of Ornstein-Uhlenbeck 
Levy processes. 
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6.1 Vectors of single and double integrals 

The following statement corresponds to Theorem 13.31 in the special case 

F={Fi,...,Fd) = {h{gi),...,h{g^),l2{hi),...,l2{K)). (32) 

The proof, which is based on a direct computation of the general bounds proved in Theorem 
13. 3| serves as a further illustration (in a simpler setting) of the techniques used throughout 
the paper. Some of its applications will be illustrated in Section 16.21 



Proposition 6.1 Fix integers n,m > 1, let d = n + m, and let C be a d x d nonnegative 
definite matrix. Let X ~ Md{^, C). Assume that the vector in h32^) is such that 

1. the function gi belongs to L^(^) n L'^{fi), for every 1 < i < m, 

2. the kernel hi € L^fJ?) (1 < i < n) is such that: (a) hi^-k^hi.^ E L^{n^), for 1 < ii,i2 1^ n, 
(b) hi € L'^{fi'^) and (c) the functions \hij^\ \hi2\ , |^iil*2l^«2l '^'^f^ I ^ii I I 1 are well 
defined and finite for every value of their arguments and for every 1 < ii,i2 ^ n, (d) 
every pair {hi, hj) verifies Assumption C, that in this case is equivalent to requiring that 



Then, 



fi{da)hf{z,a)h'^{z,a)fi{dz) < oo. 
z V Jz 



d3{F,X) < ^V^l + 52 + 53 + 54 

1 



< ^VSl + 55 + 56 + 54 
where 

m 

51 = ^ {C{ii,i2) - {gi^,gi2)L2i^^)f 

n 

52 = ^ (C(m + ji,m + j2)-2(/iji,/ij2)i2(^2))2 + 4||/ij, *^/ijJ|^2^^^ + 8||/iji^ 

ii j2=i 

m n 

■o/^/„- ™ I „•^2 , c;ll„ ,1 I, 



53 = ^j;2C(z,m + j)2 + 5||r7^*l/i,|'2 

i=i j=i 

m n 

54 = ||(7i||i3(^) +8n2^ ||/ij-||i2(^2)(||/ij||i4(^2) + V2\\hj^ /iji IIl2(;.3)) 

i=i j=i 

n 

55 = {C{m + ji,m + j2) - 2{hj^,hj^)L2(^2))^ + 4\\hj^ hj^h^ip^) x \\hj2 *i hj2\\L^(^,3) 

ii j2=i 

+ 8\\hj^ ★} /ljj|i2(^2) X ||/lj2 *1 ^^2 11^2(^2) 
m n 

Se = YY'iC{i,m + jf + 5\\g^\\l2^^)X\\hj-klhj\\L2^^2) 
i=i j=i 
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Proof. Assumptions 1 and 2 in the statement ensure that each integral appearing in the proof 
is well-defined, and that the use of Fubini arguments is justified. In view of Theorem 14.21 our 
strategy is to study the quantities in line ()15p and line ()16p separately. On the one hand, we 
know that: for 1 < i < m,l < j < n, 

DJM-)) = gi{z), -D,L-'higi{-)) = gi{z) 

DMhj{-,-)) = 2h{hj{z,-)), -D,L-'l2{hj{;-)) = h{hj{z,-)) 
Then, for any given constant a, we have: 

— for 1 < i < m,l < j < n, 

K[{a - {DJiig,,),-D,L-^higi,))f] = (a - {gi„ gi,) ^^^f ; 

— for 1 < ji, j2 < n, 

E[{a - {DMh,,), -D,L-'h{hj,))f\ 
= (a - 2(/iji,/ij2)l2(^2))^ + 4||/ij,-j *2 ^j2\\hi^J.) + ^ll^ii *i ^J2lli2(^2); 

— for 1 < i < m,l < j < n, 

E[(a - {Dj2ihj), -D,L-'h{g,))f] = a' + A\\g, A h,fm,) 

|2 



E[(a - (DJiigi), -D,L-^h{hj)) f] = + \\g, 4 h,f 



So 

1 



m = ipjsi +s2 + s^ 

where 81,82, 83 are defined as in the statement of proposition. 
On the other hand, 

Y.\D.fA = iY,\9i{z)\+2^\h{h,{z,-))\ 
i=i / \i=i j=i 

d m n 

^ \D,L-'F,\ = ^ \g,{z)\ + \h{h,{z, .))\. 
1=1 1=1 j=i 

As the following inequality holds for all positive reals a, b: 

(a + 26)2 (a + 6) < (a + 2bf < Aa^ + 326^ 
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we have, 



E 



E 



Ki=l / \i=l / 

m n \ 



i=l 



i=l 



< E 



4 +32K^|Ii(/.,(z,.))| 



m 



n 



i=i j=i 
By applying the Cauchy-Schwarz inequahty, one infers that 



^l{dz)E[\h{h{z,.))\']<JE 



Uz 



f,{dz)\hih{z,.))\^- 



X II^IIl2{/.2)- 



Notice that 



E 



f,idz)\hihiz,-m 



We have 



\i=l / \i=l / 



< 



i=l 



+8n2^||/i,-||i2(^2)(||/i,-||2,(^2) + V2||/i,4 hj\\L2 
= \\C-'\\%'\\C\\opS, 

We will now apply Lemma [2. 91 to further assess some of the summands appearing the definition 
of 52,53. Indeed, 

- for 1 < ji, j2 < n, 

ll^jl *2 ^j2\\h{fl) ^ ll^il *1 ^illlL2(;,3) X \\hj2 /ij2llL2(^3) 

— for 1 < i < m,, 1 < j < n, 

/lilli2(^) < ||5i|li2(^) X /Ij||l2(^2) 

by using the equality ||s-f^ *o 11^2(^2) = Ibf ^ Ili2(^). 
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Consequently, 

n 

S'2 < ^ {C{m + ji,m + j2) -2{hj^,hj^)L2(^^fi)f + A\\hj^-klh 
ii J2=i 



5> 
m n 



i=i j=i 
= Se 



Remark 6.2 If the matrix C is positive definite, then we have 



d2iF,X) < \\C-^\\op\\C\\l(,WSi + S2 + S3 + ^||C-i||3/2||C||op54 



-j^ 11I/2 / ^^^11 lii3/2ii 

< \\C \\op\\C\\op Si + + Sq -\ ^11^ Wop \\C\\opS4, 

by using Theorem 13.31 

The fohowing resuh can be proved by means of Proposition 16.11 or as a particular case of 
Theorem 15.71 

Corollary 6.3 Let d = m + n, with m,n > 1 two integers . Let Xq ~ A/'rf(0, C) he a centered 
d- dimensional Gaussian vector, where C = {C{s,t) : s, t = 1, . . . , d} is a d x d nonnegative 
definite matrix such that 

C{i,j + m) = 0, \/l < i < m,l < j < n. 

Assume that 

= {Fi'\ . . . ,Ff ) := {Iiig[% . . . J,{g^^))j2{hf\ . . .,l2{h^^)) 

where for all k, the kernels g^^\...,gr^ and h^\...,hn^ satisfy respectively the technical 
Conditions 1 and 2 in Proposition \6.1\ . Assume moreover that the following conditions hold 
for each k >1: 

1. 



lim E[Ff ^Ff ^] = C{s, t), l<s,t< d. 



or equivalently 



fe— !-00 

^l™^2(/i^.^\/i^.^^)i2(^2) = C(m + ji,m + j2), 1 < ji,i2 < n. 
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2. For every i = l,...,m and every j = l,...,n, one has the following conditions are 
satisfied as /c —)• oo; 

r«;ii5fiii3(,)^o;r&;ii/.fiii,(^.)^o; 

(c) \\hf 4 /if = Whf^ /if ^||l2(^3) ^ 0; 

(d) \\hfK\hf%^^.^^0. 

Then F^'^^ — )• X in law, as k — )• oo. An explicit bound on the speed of convergence in the 
distance is provided by Proposition \6.1[ 

6.2 Vector of functionals of Ornstein-Uhlenbeck processes 

In this section, we study CLTs for some functionals of Ornstein-Uhlenbeck Levy process. These 
processes have been intensively studied in recent years, and applied to various domains such 
as e.g. mathematical finance (see |15| ) and non-parametric Bayesian survival analysis (see e.g. 
[21 [16]). Our results are multi-dimensional generalizations of the content of |17^ Section 7] and 
dHl Section 4]. 

We denote by iV a centered Poisson measure over M x R, with control measure given by 
v{du), where v{-) is positive, non-atomic and cj-finite. For all positive real number A, we define 
the stationary Ornstein-Uhlenbeck Levy process with parameter A as 

= h{f^) = V2\ I I uexp(-A(t - x))N{du,dx), t > 

where f^{u,x) = -v/2Al(_oo,t](^)^6xp(— A(t — x)). We make the following technical assump- 
tions on the measure v. J^u^i'{du) < oo for j = 2,3,4,6, and f^u'^u{du) = 1, to ensure 
among other things that is well-defined. These assumptions yield in particular that 

Var(y/) = E[(y/)2] = 2X [ f u"^ exp(-2A(t - x))u{du)dx = 1 

We shall obtain Central Limit Theorems for three kind of functionals of Ornstein-Uhlenbeck 
Levy processes. In particular, each of the forthcoming examples corresponds to a "realized 
empirical moment" (in continuous time) associated with Y^, namely: Example 1 corresponds 
to an asymptotic study of the mean. Example 2 concerns second moments, whereas Example 
3 focuses on joint second moments of shifted processes. 

Observe that all kernels considered in the rest of this section automatically satisfy our 
Assumptions A, B and C. 

Example 1 (Empirical Means) 

We first recall the definition of Wasserstein distance. 

Definition 6.4 The Wasserstein distance between the laws of two W^-valued random vec- 
tors X and Y with K\\X\\j^d ,E\\Y\\^d < oo, written dw{X,Y), is given by 

d^{X,Y) = sup \E[g{X)] -E[g{Y)]\, 
where % indicates the collection of all functions g G C^(M'^) such that \\g\\Lip ^ 1- 
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1 



We define the functional A(T, A) by A(T, A) = — ^ Jq Y^dt. We recall the following limit 
theorem for A{T^\) , taken from Example 3.6 in |17) . 

Theorem 6.5 As T ^ oo, 

AiT,X) 



y/2/\ y/TT/X, 

and there exists a constant < a(A) < oo, independent of T and such that 

, (A{TA) ^\ ^ «(A) 



V V2/A 

Here, we present a multi-dimensional generalization of the above result. 
Theorem 6.6 For Ai, . . . , A^ > 0, as T ^ oo, 

A{T) = {A{T, Ai), . . . , A{T, A,)) Xb, (33) 

where Xb is a centered d- dimensional Gaussian vector with covariance matrix B = {Bij)dxd, 
with Bij = Ij ^ XiXj , 1 < i,j < d. Moreover, there exists a constant < a = a(A) = 
a(Ai, . . . , Xd) < oo, independent ofT and such that 

d,{A{T),XB)<^. 

Proof. By applying Fubini theorem on A{T, A), we have 



1 

where 



T 

Y/'dt = h[gx,T) 



2X f 

gx,T = l{-oo,T]{x)uJ ^ / exp(-A(i - x))dt 

V ^ JxVO 

E[A{T, X^)A{T,Xj)] 

I u^u{du)( I dx J exp ((Aj + Xj)x) x (l - exp(-AiT)) x (l - exp(-Ajr)) 

JM. ^J~00 Ty/XiXj 

rT 2 

+ / dx exp ((Aj + Xj)x) X ( exp(— Ajx) — exp(— AjT)) x ( exp(— A^x) — exp(— AjT) 

Jo T^yXiXj 

^ X (l - exp(-Air)) X (1 - exp(-Ajr)) + T - ^ x (1 - exp(-Air)) 



Ty'^AiAj ^Aj + Xj Xi 
-^(1 - exp(-A,r)) + Y^i^ - exp(-(A. + A,)r) 

Aj Ai -\- Aj 

2 /l\ 

: + - as OO. 



y/XiXj I r 
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And we may verify that ||5A.T||^3(^j^f;2,) ~ y^' A G R. (See [ITj and [l8] for details.) 

Finally, we deduce the conclusion by using Corollary 14.31 ■ 



Example 2 (Empirical second moments) 

We are interested in the quadratic functional Q{T, A) given by: 

Q(r, A) := Vf (^^ [o't^)^'^^ - 1) ' r > 0, A > 

In |17) and [18] . the authors have proved the following limit theorem for Q(T, A). (See Theorem 
7.1 in [17] and Proposition 7 in [18J) 

Theorem 6.7 For every A > 0, as T ^ oo, 



Q{T,\) -.= ^{1 j\Y,^fdt-l^ 



where X ~ A/'(0, 1) is a standard Gaussian random variable and = u'^i^^du) is a constant. 
And there exists a constant < /3(A) < oo, independent ofT and such that 

We introduce here a multi-dimensional generalization of the above result. 
Theorem 6.8 Given an integer d> 2. For Xi, . . . , X^i > 0, as T ^ oo, 

Q{T) = iQ{T, Ai), . . . , Q{T, A,)) Xc, (34) 

where Xq is a centered d- dimensional Gaussian vector with covariance matrix C = {Cij)dxd> 
defined by 

4 



Cij = - — — + c^, l<i,j <d, 



Aj + Aj 

and cl = JjgU^i^(du). And there exists a constant < /3(A) = /3(Ai, . . . , A^) < oo, independent 
of T and such that 

/3(A) 

dsiQiT),Xc) < 



/T 

Proof. For every T > and A > 0, we introduce the notations 

/ / / -'-(-oo T)2 (2^) 2;') / , / \ / \ 

Hx,t{u,x;u ,x) = (uxu) exp (A(x + x )) x (1 - exp(-2Ar)) x l(^va:'<o) 



+ exp (A(x + x')) X ( exp(-2A(x V x')) - exp{-2XT)) x l(^va: 



'>o) 



HIt{u,x) = ^2 l( °y)^^) (^exp(2A:r) x (l - exp(-2Ar)) x l(,,<o) 



+ exp(2Ax) X (exp(-2Ax) - exp(-2Ar)) x 1(, 



x>0) 
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By applying the multiplication formula (|6]) and a Fubini argument, we deduce that 

Q(T, A) = hiVTHlT) + I2{VtHx,t), 

which is the sum of a single and a double Wiener-Ito integral. Instead of deducing the 
convergence for {Q{T, Ai), . . . , Q{T, A^)), we prove the stronger result: 

ihiVfHl^T), hiVTHl^j,), UVtHx^t), l2iVfHx,,T)) Xd (35) 

as T — )■ oo. Here, is a centered 2d-dimensional Gaussian vector with covariance matrix D 
defined as: 



' cl, ifl<i,j<d 



-, if d+l <i,j <2d 



Aj + Aj 
^ 0, otherwise. 

We prove ()35p in two steps (by using Corollarv l6.3p . Firstly, we aim at verifying 



lim E[fPfP] = D{i,j), l<i,j < 2d, 



r-s>oo 



for 



F, 



(T) 



\ hiVTHx^r), if d + 1 < A; < 2d 
Indeed, by standard calculations, we have 



T / H^,rp{u,x)H^.rp(u,x)i'{du)dx 



1 



1 



T n2(Ai + A,: 

-X h - 

2A,- 



;i - exp(-2Air)) X (1 - exp(-2Ajr)) + T - — x (l - exp(-2Air)) 

2 A,; 



X (1 - exp(-2A,r)) + ^p^-^l^ X (1 - exp(-2(A, + A,)r) 



as r — )• oo, 



and 



2T / Hx-^t{u,x;u' ,x')Hx.t{u,x;u' ,x')i'{du)i'{du)dxdx' 



2^( l-exp(-2AiT))x(l-exp(-2A,r)) , 2 1 , n\ 



-^x(l-e,p(-2A,r))+^^^ 



(1 - exp(-2(A, + Aj)T) 



Secondly, we use the fact that for A = Ai, . . . , A^, the following asymptotic relations holds as 
oo: 
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(a) llV^^A,rlli3(di.da;) ~ 



1 

1 

1 

Hx,t) *i (VTi7A,T)||L2((rfj,rf2.)3) ~ - 

1 



(c) IKVtFa.t) 4 {VTHx,T)\\LHdud.) = \\iVTHx,T) A iVTHx,T)\\LHidud.)^) ~ 

(d) ||(\/rFA,T) *l iVTHx,T)\\LH(dud.r) ~ 



(e) \\{VTHy)4{VTHx,T)\\ 



/T 

The reader is referred to [17| Section 7] and [18i Section 4] for a proof of the above asymptotic 
relations. ■ 

Example 3 (Empirical joint moments of shifted processes) 

We are now able to study a generalization of Example 2. We define 

Qh{T, A) := Vt (^i £ Yt%li,dt - exp(-A/i)^ , /i > 0, T > 0, A > 0. 

The theorem below is a multi-dimensional CLT for Qh{T,X). 
Theorem 6.9 For Xi, . . . , > and h > 0, as T ^ 00, 

Qh{T) = {Qh{T, Ai), . . . , Qh{T, Xd)) Xe, (36) 

where Xe is a centered d- dimensional Gaussian vector with covariance matrix E = {Eij),^^^^, 
with 

Eij = ^ _^ ^ + cl exp {Xi + Xj)h^ , l<i,j<d 

and = J^u'^i'{du). Moreover, there exists a constant < 7(/i, A) = 7(/i, Ai, . . . , Xd) < 00, 
independent of T and such that 

ds{Qh{nxE)<^^^ 

Proof. We have 

fv^^Y^X^dt = f \^{f^)h{f^_^^)dt 
Jo Jo 

= £ ' {hift' 4 ft\h) + hUt A ft+h) + ft A ft\h)dt 

= £ \h{hl^) + h{h;i^) + c)dt 

= h{THlT) + h{TH*^^) + C{T-h) 

by using multiplication formula (j6|) and Fubini theorem. By simple calculations, we obtain 
that 

'ht,h(.u,x;u',x') = 2Al(_oo,t]x(-oo,t+fc](2;,2;') x W exp(-A(2t + h-x- x')) 
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kl'^{u,x) = 2Al(_oo,t](2;) xu'^exp{-X{2t + h-2x)) 



as well as 



h*^^{u, x)dt 



w 



i(-oo,T](2;) 



T 

+1|,<0} X (1 - exp(-2Ar))) 



X exp(A(2x — h)) x M|^>o} ^ (exp(— 2Ax) — exp(— 2AT)) 



HxAu,x]u',x 



1 



h^^{u, x; u ,x')dt 



^ ^^^^^^^ ^ ^, _ 



X {l{xv{x'~h)>o} X (exp(-2A(x V {x - h))) - exp(-2Ar)) 
+l{xv(x'-h)<o} X (1 - exp(-2Ar))^ 
Similar to the procedures in the precedent example, we prove the stronger result: 

ihiVTHlf^^), . . .,h{VTHlf^^),h{VTH^,^^T), . . .,I2{Vth',^^t)) X^. (37) 

as T — >• oo. Here, Xj^h is a centered 2d-dimensional Gaussian vector with covariance matrix 
D'' defined as: 



^ clexp{-{Xi + Xj)h), iil <i,j<d 
4 

ifd+l<i,j <2d 
otherwise. 



Aj + Aj 



0, 



We have 
T 



H*^^{u, x)H*^^{u, x)v{du)dx 



1 / /-o 

= -cin dxexp ((Ai + Aj)(2x - h)) x (l - exp(-2Air)) x (l - exp(-2Ajr)) 

+ / (ij;exp ((Aj + Aj)(2x - h)) x (exp(-2Aij;) - exp(-2Air)) x (exp(-2Ajx) - exp(-2Ajr) 
Jo 



exp(-(Ai + \j)h) + ^ ( ^ 1 • as r — )> oo. 



We notice that 



H^rpiu, x; u , x) = Hx^t{u, x; u ,x' — h) 



Then, as shown in the proof of Theorem 16. 8| we have 

2T [ HIt{u, x)hIt{u, x)u{du)dx = -4l- + o(\; 



as T — >■ oo. 
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(6) WVfHl^Wl^^^^^^^^^^^—; 



Just as the precedent example, we may verify that for A = Ai, . . . , and h > 0, the following 
asymptotic relations holds as T — )• oo: 

(a) WVfHlpll,^^^^^^ ~ 



(c) WiVTHl^) 4 {VTHI^\\l2^,,,,) = \\iVTHx,T) 4 (v/Ti/^^_^)|U2((,,,,)3) ~ 

(e) ||(x/r/?*;^) {Vf''H,^T)\\LHaud.) ~ 

We conclude the proof by analogous arguments as in the proof of ([3; 



The calculations above enable us to derive immediately the following new one-dimensional 
result, which is a direct generalization of Theorem 5.1 in |17) . 

Corollary 6.10 For every A > 0, as T ^ oo, 

Qh{T, A) ^^+c2exp(-2A/i) x X 

where X ~ A/'(0, 1) is a standard Gaussian random variable. Moreover, there exists a constant 
< 7(/i, A) < oo, independent ofT and such that 

dJ^^^^=,x]<l^ 

\^V2/A + c2exp(-2A/i)' J' VT 
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7 Appendix: Malliavin operators on the Poisson space 



We now define some Malliavin-type operators associated with a Poisson measure N, on the 
Borel space {Z,Z), with non-atomic control measure /x. We fohow the work by Nualart and 
Vives |14| , which is in turn based on the classic definition of Malliavin operators on the Gaus- 
sian space (see e.g. [7t 

(I) The derivative operator D. 

For every F G L^(F), the derivative of F, DF is defined as an element of L2(P;L2(^)), 
that is, of the space of the jointly measurable random functions u : x Z i— )■ R such that 
<oo. 

Definition 7.1 1. The domain of the derivative operator D, written domZ?, is the set of 
all random variables F S L'^(P) admitting a chaotic decomposition ([ip such that 

^kk\\\fk\\l2i^^k) < oo, 

k>l 

2. For any F G domD, the random function z i— >• D^F is defined by 

oo 

D,F = Y,kIk-i(.fk(.z,-)). 

k>l 

(II) The divergence operator 5. 

Thanks to the chaotic representation property of TV, every random function u £ L^(P, L'^{fi)) 
admits a unique representation of the type 

oo 

u, = Y,Wk{z,-)), zeZ, (38) 

fc>0 

where the kernel is a function of A: + 1 variables, and fk{z, •) is an element of L1{^^). The 
divergence operator 5{u) maps a random function u in its domain to an element of L^(P). 

Definition 7.2 1. The domain of the divergence operator, denoted by dova.6, is the col- 
lection of all u G L^(P, L^(/i)) having the above chaotic expansion h3^) satisfied the 
condition: 

E(^+i)'ii/fciii2{/.(fc+i))<°°- 

2. For u G dom(5, the random variable 5{u) is given by 

^(u) = ^h+iifk), 
k>0 

where f^ is the canonical symmetrization of the k + 1 variables function fi^. . 
As made clear in the following statement, the operator 6 is indeed the adjoint operator of D. 
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Lemma 7.3 (Integration by parts) For every G G domZ) and u G dom(5, one has that 

nG6{u)]=n{DG,u)L^^)\. 
The proof of Lemma 17.31 is detailed e.g. in |14) . 

(III) The Ornstein-Uhlenbeck generator L. 

Definition 7.4 1. The domain of the Ornstein-Uhlenbeck generator, denoted by domL, is 
the collection of all F G L^(P) whose chaotic representation verifies the condition: 

fc>l 

2. The Ornstein- Uhlenbeck generator L acts on random variable F G domL as follows: 

LF = - Y^khifk). 

k>l 

(IV) The pseudo-inverse of L. 

Definition 7.5 1. The domain of the pseudo-inverse of the Ornstein-Uhlenbeck generator, 
denoted by L^^ , is the space Lq{F) of centered random variables in L^(P). 

2. ForF=Y. hifk) G Ll{¥) , we set 

k>l 

L~'F = -Y,\h{f,). 

k>l 
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